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We study the interaction of a weak probe field, having two orthogonally polarized components, 
with an optically dense medium of four-level atoms in a tripod configuration. In the presence 
of a coherent driving laser, electromagnetically induced transparency is attained in the medium, 
dramatically enhancing its linear as well as nonlinear dispersion while simultaneously suppressing 
the probe field absorption. We present the semiclassical and fully quantum analysis of the system. 
We propose an experimentally feasible setup that can induce large Faraday rotation of the probe 
field polarization and therefore be used for ultra-sensitive optical magnetometry. We then study the 
Kerr nonlinear coupling between the two components of the probe, demonstrating a novel regime of 
symmetric, extremely efficient cross-phase modulation, capable of fully entangling two single-photon 
pulses. This scheme may thus pave the way to photon-based quantum information applications, such 
as deterministic all-optical quantum computation, dense coding and teleportation. 

PACS numbers: 42.50.Gy, 07.55. Ge, 03.67.-a 



I. INTRODUCTION 

Electromagnetically induced transparency (EIT) in 
atomic media is a quantum interference effect that results 
in a dramatic reduction of the group velocity of propa- 

Stingprobe field accompanied by vanishing absorption 
the quantum interference is usually very 
sensitive to the system parameters, various schemes ex- 
hibiting EIT are attracting growing attention in view of 
their potential for significantly enhancing nonlinear op- 
tical effects. Some of the most representative examples 
include slow-light enhancement of acusto-optical interac- 
tions in doped fibers , trapping light in optically dense 
atomic and doped solid-state media by coherently con- 
verting photonic excitation into spin excitation 's','^,'?! or 
by creating a photonic band gap via periodic modulation 
of the EIT resonance |^, and nonlinear photoi>photon 
coupling using N configuration of atomic levels ^ [13 • 

EIT is based on the phenomenon of coherent popula- 
tion trapping 1 , 2] , in which the application of two laser 
fields to a three-level A system creates the so-called "dark 
state" , which is stable against absorption of both fields. 
Dark states are also found in several other multilevel sys- 
tems, one of them being four-level atoms interacting with 
three laser fields in tripod configuration. Tripod atoms 
proved to be robust systems for "engineering" arbitrary 
coherent superpositions of atomic states using an ex- 
tension of the well-known technique of stimulated Raman 
adiabatic passage (STIRAP) 12]. Parametric genera- 
tion of light in a medium of tripod atoms, prepared in 
a certain coherent superposition of ground states, has 
been recently discussed in . In a related work, it was 
shown that enhanced nonlinear conversion between two 
laser pulses is attainable in a medium of A atoms with 
spatially dependent ground state coherence 0|. In the 



present paper we undertake a detailed study of propaga- 
tion of a weak probe field through a medium of tripod 
atoms under the conditions of EIT 0|. We show that 
this system can support large magneto-optical rotation 
(MOR) of the probe field polarization, accompanied by 
negligible absorption. It can therefore be used for ultra- 
sensitive optical magnetometry, with the sensitivity com- 
parable to (or better than) other hitherto studied MOR 
schemes In contrast to these schemes, where the 

basic mechanism of nonlinear MOR is the probe field in- 
duced coherence between the Zeeman sublevels of atomic 
ground state 0, 0| , in our case the MOR results from 
an extraordinary dispersion induced by a strong driving 
field in the EIT regime. Hence, by simply changing the 
intensity of the driving field, one could control the po- 
larization rotation of the weak probe field. We note that 
an interferometric measurement of the magnetic field in- 
duced phase shift of the probe, subject to EIT in the 
presence of a driving field, can yield sensitivity of the or- 
der of lO"'^^ G These studies and our present con- 
tribution reveal the significant potential for improving 
the sensitivity of Faraday magnetometers to small mag- 
netic fields as compared to conventional optical pumping 
magnetometers 20| . 

Another motivation for the present work is its rele- 
vance to the field of quantum information (QI), which 
is attracting broad interest in view of its fundamental 
nature and its potentially revolutionary ap plic ations to 
cryptography, teleportation and computing |2l|. Among 
the various QI processing schemes of current interest 
mm mill Eg, those based on photons lUmi have 
the advantage of using very robust and versatile carriers 
of QI. Yet the main impediment towards their opera- 
tion at the few-photon level is the weakness of optical 
nonlinearities in conventional media As mentioned 
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FIG. 1: Level scheme of tripod atoms interacting with weak 
probe £ and strong driving Ed fields. Upper inset: coprop- 
agating probe with circularly left- and right-polarized com- 
ponents £1,2, and driving Ed fields pass through the atomic 
medium that is subject to the longitudinal magnetic field B. 
Lower inset: Perpendicular arrangement of the probe and 
driving fields, that is suitable for cold atomic gas. 
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above, EIT schemes with atoms having N configuration 
of levels have opened up a possibility of achieving en- 
hanced nonlinear coupling of weak quantum fields at the 
single-photon level [9||, JO]. The main hindrance of such 
schemes is the mismatch between the group velocities of 
the pulse subject to EIT and its nearly-free propagat- 
ing partner, which severely limits their effective interac- 
tion length This drawback may be remedied by us- 
ing an equal mixture of two isotopic species, interacting 
with two driving fields and an appropriate magnetic field, 
which would render the group velocities of the two pulses 
equal [25. Here we propose an alternative, simple and 
robust approach which relies solely on an intra-atomic 
process, without resorting to two isotopic species and us- 
ing just one driving field [Tsl l29l |. In our scheme, two 
orthogonally polarized weak (quantum) fields, acting on 
adjacent transitions of tripod atoms, propagate with the 
same group velocity and impress large conditional phase 
shift upon each other. 

The paper is organized as follows. In Sec. ^1 we for- 
mulate the theory and give an analytical solution of the 
equations of motion for the two components of the weak 
probe field. In Sec. Illll we discuss the setup and sen- 
sitivity limits of the optical magnetometer. Section IIVI 
is devoted to the study of feasibility of strong nonlinear 
interaction and entanglement between two orthogonally 
polarized weak quantum fields, aimed at quantum infor- 
mation applications. Our conclusions are summarized in 
Sec. El 

II. FORMULATION 

We consider a near-resonant interaction of two opti- 
cal fields with a medium of atoms with tripod configu- 
ration of levels (Fig. QJ. The medium is subject to a 
longitudinal magnetic field B that removes the degener- 
acy of the ground state sublevels. The Zeeman shift of 
levels |1) and |2) is given by HA = fiBMpg^B, where 
fj,B is the Bohr magneton, (/f is the gyromagnetic factor 



FIG. 2: Absorption and dispersion spectra (S = LOp — W31) for 
the £1 and £2 components of the probe in the presence of a 
strong driving {fid = O.BF) and a weak magnetic (A = O.IF) 
fields, in units of the linear resonant absorption coefficient ao. 



and Mp = ±1 is the magnetic quantum number of the 
corresponding state. All the atoms are assumed to be 
optically pumped to the states |1) and |2) which thus 
have the same incoherent populations equal to 1/2. A 
linearly polarized weak (quantum) probe field £ has a 
carrier frequency ujp and wavcvector kp parallel to the 
magnetic field direction. Its two circularly left- and 
right-polarized components £1^2 ^ct on the atomic tran- 
sitions |1) 1 3) and |2) |3), with the detunings 
Si,2 = u>p—u>2i—kpV^A, where uj^i = 10^2 is the frequency 
of the unshifted atomic resonance and kpV is the Doppler 
shift for the atoms having velocity v along the probe field 
propagation direction. A strong classical cw field Ed, 
having frequency Ud and wavevector kd — kp, is driving 
the atomic transition |3) <-> |4) with the Rabi frequency 
— paiEd/h, where p^i, is the dipole matrix element 
on the transition \fj,) In the coUinear Doppler- 

free geometry shown in Fig. ^ upper inset, the driving 
field has to be circularly left or right polarized, in order to 
couple to a single magnetic sublevel |4). Its Zeeman shift 
HA' — fiBMprgpiB is incorporated in the detuning of the 
driving field via 6d — ujd — — kdV + A', where 0^34 is 
the atomic resonance frequency for zero magnetic field. 
Note that in the case of cold atomic sample (Doppler 
broadening of the atomic resonance is smaller than the 
ground-state spin relaxation rate), one can employ the 
perpendicular geometry of Fig. ^ lower inset, where the 
driving field is linearly tt polarized while the Zeeman shift 
of level |4) vanishes. A' = since Mpi = 0. 

To illustrate the scheme, we plot in Fig. |21the absorp- 
tion and dispersion spectra of the two components of the 
probe field £ for the case Sd = 0. In the presence of mag- 
netic field, the spectra for the £1 and £2 are shifted with 
respect to each other by the amount equal to the Zeeman 
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shift 2A. When the probe field is resonant with the un- 
shifted (A = 0) atomic transitions, Wp = uj^^ = Wgj, due 
to the steep and approximately linear slope of the disper- 
sion in the vicinity of (5i,2 = 0, upon propagating through 
the medium the two components of the probe experience 
equal and opposite phase shifts (pi = — 02 which results in 
a net polarization rotation of the field, $ = \{4>2 — <t'i)- 
If the Zeeman shift is small compared to the width of 
the EIT window for both components of the probe, the 
absorption remains much smaller than the phase shift. 
Thus, a weak magnetic field can induce an appreciable 
polarization rotation accompanied by negligible absorp- 
tion, allowing for extremely sensitive magnetometry (Sec. 
mill . In addition to the large linear phase shift, each com- 
ponent experiences a nonlinear cross-phase modulation. 
Although this cross-phase modulation is typically small 
compared to the linear phase modulation, it is neverthe- 
less several orders of magnitude larger than that in con- 
ventional media 9] . It can therefore be used for quantum 
information applications based on photon-photon inter- 
action and entanglement (Sec. IIV|I . 

Let us now consider the scheme more quantitatively. 
We describe the medium using collective slowly vary- 
ing atomic operators (j^^{z, t) = ^ \ ' a-ver- 
aged over small but macroscopic volume containing many 
atoms A^^ — {N/L)dz ^ 1 around position z, where TV 
is the total number of atoms and L is the length of the 
medium . The two components of the quantum probe 



field are described by the corresponding field operators 
£i.2- In a frame rotating with the probe and driving field 
frequencies, the interaction Hamiltonian has the follow- 
ing form 

N 

H = h— I dz[Si&ii + 62(722 + Sd^44 

^ Jo 

-.9(^10-31 +^20-32) - ^d^34 +H.C.]. (1) 



Here g = psi y/ojp / {2heoAL) , with A being the cross- 
sectional area of the probe field, is the atom-field coupling 
constant, which is the same for both circular components 
£1^2 due to the symmetry of the system (|p3i| = |p32| 
while the opposite signs of the Clebsch-Gordan coeffi- 
cients on the transitions |1) |3) and |2) |3) can 
be incorporated into the atomic eigenstate via the trans- 
formation |1) Using the slowly varying enve- 
lope approximation, we obtain the following propagation 
equations for the quantum field operators 

(^■^^+c^^£iiz,t) = igNaia, (2a) 

f^+c^j^2(^,t) = ^gNa23■ (2b) 

The equations for the atomic coherences are given by 



d_ 

di 
d_ 
dt 

d . 

d „ 
d . 



fl3 



i{6i - 62) - 7c]'5'i2 - ig£i^32 + ig£l^i3 + A2, 

- j <5-i3 + ig£i{^ii - 0-33) + ig£2(yi2 + ifld^u + F13, 
iiSi - 6d) - 7c]cti4 - ig£id-34 + i^*d^i3 + A4, 

«<52 - - j (723 + W^2{022 - <5'33) + W^1^21 + i^d^24 + F23, 
[i{62 - 6d) ~ lc]cf24 - ig£2Cf34 + i^*d^23 + F24, 

- (iSd + ^ j '5-34 - ■ig£l^i4 - ig£l^24 + i^*di^33 - 0-44) + ^4, 



(3a) 
(3b) 
(3c) 
(3d) 
(3e) 
(3f) 



r 



where 7^ is the ground state coherence (spin) relaxation 
rate, F is the decay rate of the excited state |3) and 
Ffj_i, are (5-correlated noise operators associated with the 
relaxation. 

We now outline the solution of Eqs. © in the weak 
field limit. To this end, we assume that the Rabi fre- 
quencies g£i^2 of the quantum fields are much smaller 
than ^^d and the number of photons in £12 is much less 
than the number of atoms, therefore an = (T22 — 1/2 
while a'33 = 1744 = (734 ~ 0. We may thus treat the 



atomic equations perturbatively in the small parameters 
g£i,2/^d- In the first order, from l)3b|l and (j3d|l we have 



g£i 
'2nd' 



^(1 



9£2 

2Qd 



Substituting these into l|Sc|) and l|Se|l . and neglecting for 



4 



now the spin relaxation, we obtain 



'13 



^(1) 



d_ 

dt 



ig£i ggi((5i - Sd) 

2\nd\^ 

5^2 ((52 - Sd) 



2\nd\^ 



In these equations, the last equalities result from the 
adiabatic approximation, i.e., we assume that the probe 
pulse changes slowly enough so that the atoms follow the 
field adiabatically. Quantitatively, the adiabatic evolu- 
tion requires that the rate of change of the probe field, 
max[9t£/£] ~ T~^, where Tp is the temporal width of 



the pulse, should be smaller than any transition rate be- 
tween the system's eigenstates, so that no nonresonant 
transition is induced 0, I^Ol ■ 

We next write Eq. (|3a|l in an integral form and perform 
the integration. 



CTl2 



.9ticr^2 



9^2^ 



t;^(l) 



2"13 



»7c 



2A 



^— i(2A-i7e)t 



where a ~ [1 + (TpA)-^]^^ is the adiabaticity parameter. 
Thus in the adiabatic limit Tp ^ as well as for 

times t ^ (for any A), the term proportional to a 
vanishes. Substituting the above expressions into 



(7i4 = 


5^1 

2nd 


9^2 . 


i 

n'd 




^^ 




~ Fi3 


(724 = 


g£2 
2nd 


9^1 - 


i 

n'd 




^^ 

" 2y 


41' 


— F23 



after some algebra, we finally arrive at the following set of equations 



an 



0'23 



0'24 



n* 

9^1 
2nd 

i 

n* 

9^2 
2nd 



<T14 + ^i^l4, 



iSi+iT/2)iS,-6d) , 9^122 A 



\nd\^ii-/c-2A) 



dt 



i{S2 - Sd) + 7c 



0'24 + TT^F24, 



{S2+iT/2)iS2~Sd) g'/i2A 

|f^dP(z7c + 2A) 



?r^23, 



(4a) 
(4b) 
(4c) 
(4d) 



where Ij = SjSj is the dimensionless intensity (photon-number) operator for the jth field. 

From now on we focus on the case of usp — ui^^ = uj^2- Substituting Eqs. Q into Eqs. the equations of motion 
for quantum fields are obtained as 



d 


1 


d 


dz 


^ — 


di 


d 


1 


d 


dz 


Vn 


dt 



ii = -mil - i{A + Ad){si - vil2)£i + A, 

£2 = ~K2£2+i{A- Ad){s2-V2h)£2+:F2, 



(5a) 
(5b) 



where A^ = ujd — + A' is the driving field detuning, lation coefficients, 



Kl,2 



S1.2 



2c\nd\^ 
Ng' 



7c + 
1 



r(A±Ad)2 



A(A±Ad) 



^71,2 



Ng^2A 



2c|^!d|4(2AT^7c) 



are, respectively, the linear absorption and phase modu- 



are the cross-coupling coefficients, Vg^'^^ = (1/c-f 51,2) ^ 
are the group velocities of the corresponding fields, and 
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are the noise operators having the properties 

{P^{z)T]{z')) = 2k,%<5(z-z'). 

In deriving Eqs. (O, we have assumed that the usual 
EIT conditions \Q.d\'^ > (A ± Ad)A:p^d'u, 7c(r + kp^^v), 
where v is the mean thermal atomic velocity, are satisfied, 
allowing us to neglect the Doppler induced absorption. 
On the other hand, since the terms containing kpV enter 
Eqs. linearly, the net phase-shift of the quantum 
fields, due to the Doppler shifts of the atomic resonance 
frequencies, averages to zero. Note also that if states 
|2) and |4) belong to different hyperfine components of a 



common ground state, the frequencies ujp and lo^ of the 
optical fields differ from each other by at most a few GHz, 
LOp — LOd — uj^i ^ ^^p,d- Then, as seen from Eqs. Q, the 
difference (fcp — kd)v in the Doppler shifts of the atomic 
resonances |2) |3) and |4) |3) is negligible. 

When A(A ± A^) <C l^^dP, the group velocities of Ei 
are £2 are practically the same, Wg^'^"* — Vg. Expressing 
the atom-field coupling constant g through the linear res- 
onant absorption coefficient ap = p\^ujpp/ (hctoT) for the 
transitions |1), |2) |3) as Ng^ = a^cT /2 and assuming 
that the density of atoms p — N/{AL) is large enough 
so that aocT > 4|r2dp, we have Vg ~ 4:\fld\'^ / (aoT) < c. 
Then the solution of Eqs. can be expressed in terms 
of the retarded time t = t — z/vg as 



£i{z,t) ~ £i(0,r)exp —Kiz + i(f)i{z,0,t) + / dz'J^i{z') exp ~ki{z — z') + i(f>i{z, z' ,t) 
L J Jo L 

£2iz,t) — £2(0,r)exp —K2Z + i4>2iz,0,t) + / dz'T2{z')exjp —K2{z — z') + i(j>2{z,z',t) 
L J Jo L 

where the phase operators are given by 

4>i(z,z',t) = -si(A + Arf)(z-z') + r;i(A + A<j) f dz" i2{z" ,t + z" /vg), 
^2{z,z',t) = S2(A- Arf)(z-z')-r72(A- Ad) / dz" h{z" ,t + z" /vg). 

J z' 

I 



(6a) 
(6b) 



These are the central equations of this paper. The first 
terms in Eqs. JHl describe the linear attenuation and the 
phase shift of the corresponding quantum field E\ 2 upon 
propagating through the medium, while the second terms 
account for the noise contribution. Note that although 
the expectation values of the field operators decay, yet 
slowly, with the propagation, due to the presence of the 
noise operators, their commutators are preserved 0|- We 
emphasize again that Eqs. © are obtained within the 
weak field and adiabatic approximations. 

In the following section we explore the classical limit 
of Eqs. © for the purpose of sensitive magnetometry. 
In Sec. IIVI we study the quantum dynamics of the sys- 
tem and show that our scheme is capable of realizing 
strong nonlinear interaction and entanglement between 
two tightly focused quantum fields at the single-photon 
level. 



III. OPTICAL MAGNETOMETER 

Let us consider the classical limit of Eqs. l|njl, by 
replacing the operators E\,2 with the corresponding c- 
numbers E\,2 and dropping the noise terms. The equa- 
tions for the two circularly polarized components of the 



cw probe field have the form 

£i(z) = £i(0)e-'^i^e^'^i(^\ (7a) 

£2[z) = 52(0)e-"^"e^'^^("\ (7b) 

where the absorption coefficients and phase shifts can be 
expressed through the group velocity Vg as 

r(A±Arf)2 



7c 

Kl,2 = — 



Vg\^d\' 



A. 



A(A 



A + Arf 52/2(0) 1 



A- A„ 



W 2fc2 

A(A-A,)2 



Vg Vg\nd\^ 

A- Ad 5^/1(0) 1 -e-2fci^ 



(8a) 



(8b) 



(8c) 



Vg \nd\^ 2fci 

When the absorption is small, ki,2^ <C 1, z G {0,L}, 
which requires that Vg/^c > L and A^-l- A^ < jc\^d\'^/^, 
the polarization rotation of the probe field <&(z) — 
5 [02 (2) - </'i(^)] is given by 

A .g2/(o) 



A A(A2 

— z+ ^ 



A?) 



Vgl^dl'^ 



Vg \nd\^ 



(9) 
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where 7(0) = /i(0) = /2(0) since the probe is Hnearly 
polarized at the entrance to the medium. In Eq. Q , the 
first term is hnear in the magnetic field while the second 
term has a cubic dependence on the field strength. Here 
we focus our attention on the measurement of dc mag- 
netic fields employing the dominating linear term. We 
wish, however, to point out that the presence of the cubic 
term may facilitate the detection of ac fields oscillating 
slowly compared to the bandwidth of the magnetometer, 
which is limited by the bandwidth of the EIT window 

m 



r y/SiTpL 



(10) 



Then the spectrum of along with the fundamental fre- 
quency of the magnetic field, will also contain its third 
harmonic which, for very small frequencies, may be eas- 
ier to detect jS^I- This issue is beyond the scope of this 
paper and will be addressed elsewhere. Finally, the last 
term of Eq. ((SJ, being proportional to the product of 
the magnetic field strength and probe intensity, is a con- 
sequence of Kerr-type nonlinear interaction between £i 
and £2, which is the subject of the following section. 

We consider a magnetometer setup in the "balanced 
polarimeter" arrangement [l^ . in which, at the exit from 
the medium z = L, a polarizing beam splitter oriented 
at 7r/4 to the input polarizer [$(0) = 0] is used as an 
analyzer. Then the detector signal S is represented by 
the difference of photocounts in the two channels of the 
analyzer 



S = 2n\r, 



-2kL 



sin[$(i)]cos[$(i)]. 



(11) 



where riin — PijJ,ml {^p) — '21{0)ctm/ L, with being 
the input power of the probe, is the number of photons 
passing through the medium during the measurement 
time tm- For simplicity, we neglect the difference be- 
tween the absorption coefficients for the circularly left- 
and right-polarized components, ki ~ K2 = k, which 
amounts to neglecting the cUipticity of the output field 
(e ~ 1) since = 2rAAdL/{vg\nd\^) < 1. 

The most important characteristic of a magnetometer 
is its sensitivity to weak magnetic fields, which is lim- 
ited by the measurement noise. The smallest detectable 
magnetic field Bmin can be defined as being the field for 
which the signal is equal to the noise. In our system, 
the total noise M = Mat + Mhot has two contributions, 
atomic noise A/'at and photon counting shot-noise A/'shot ■ 
The atomic contribution is due to the spontaneous pho- 
tons reaching the detector during the measurement time. 



A/'at = r(a33)iVi„ 



A 



47rL2' 



where the detector area is assumed to be equal to A. 
For vanishing magnetic field A < 7c, we have (o'33)iV ~ 
aocr7c/(0)/(4|rirf|^) and the atomic noise is given by 



327r|f7<7|4L' 



For physically realistic parameters (see below), the 
atomic noise term is small compared to the photon count- 
ing shot-noise 



shot 




In the limit of weak magnetic field, retaining only the 
linear in magnetic field term, from S ~ 2riin$ > A4hot 
we obtain 



> 



(12) 



For realistic experimental parameters, Wp = 3 x 10^^ 



rad/s, r = 10^ s- 



P 



10 



13 



(ao ~ 10^ cm^i). 



|5f| = 1/2, Od ~ r, L = 10 cm, Pjn = 1 mW, <,„ = 1 s, 
the minimum detectable magnetic field Bniin 5. 10"^^ G, 
which is of the same order as that of Thus, 
concerning the magnetometer sensitivity, our scheme is 
essentially equivalent to the one proposed in |19|| . where 
an interferometric measurement of the magnetic field in- 
duced phase shift of a probe field, subject to EIT with 
A atoms, was studied. Experimentally, however, measur- 
ing the polarization rotation of the probe, as suggested 
here, may be more practical than measuring its phase 
shift in the setup of J^] , which employs a Mach-Zehnder 
interferometer. 



IV. CROSS-PHASE MODULATION 

In order to rigorously describe the nonlinear interac- 
tion between the weak pulsed fields, we now turn to the 
fully quantum treatment of the system. When absorp- 
tion is small enough to be neglected, from Eqs. © we 
have 

£i{z,t) = fi(0,T)exp[ir;(A + Arf)f|(0,T)52(0,T)z], 

(13a) 

£2{z,t) = fi(0,T)expH7?(A- A<i)£|(0,r)fi(0,r)z], 

(13b) 

where the cross-phase modulation coefficient is given 
by ^7 = 5^/(^ffl^<iP) (assuming ^ A), while the 
linear phase-modulation is incorporated into the field 
operators via the unitary transformations £i{z,t) — > 
£i(z,i)e"*^i('^+'^<')^ and 4(^,i) ^ 4(2, i)e'^"''^~^'''^ 
These traveling-wave electric fields can be expressed 
through single mode operators as £j{z, t) = aj(t)e*''^ 
(j ~ 1,2), where is the annihilation operator for the 
field mode with the wavevector hp + q. The single-mode 
operators a"^ and a'^ possess the standard bosonic com- 



mutation relations [aj 



The continuum 



of modes scanned by g G {—5q/2,5q/2} is bounded by 
the EIT window via 6q < Suj/c j2^. The finite quanti- 
zation bandwidth Sq for the field operators leads to the 
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equal-time commutation relations 

= %^sinc[<5q(z- z')/2], 

where sinc(a:) — sin(a;)/x. 

Before proceeding, we note that Eqs. Hl^lla re simi- 
lar to the corresponding equations of Ref. [23|, where 
the cross-phase modulation between two quantum fields 
was mediated by atoms with N configuration of levels 
, while the group velocity mismatch between the fields 
was compensated by using a second kind of A-atoms con- 
trolled by an additional driving field. In contrast, our 
scheme relies solely on an intra-atomic process employ- 
ing only one driving field that causes simultaneous EIT 
for both fields and their cross-coupling. It is therefore de- 
prived of complications associated with using mixtures of 
two isotopic species of atoms or invoking cavity QED 
techniques [2y|. 

The most classical of all the quantum states is the 
coherent state. To compare the classical and quantum 
pictures, we therefore consider first the evolution of in- 
put wavepacket jV'in) = l^i) ^ \a2) composed of the 
multimode coherent states \aj) = Ilg jap (j — 1,2). 
The states \aj) are the eigenstates of the input oper- 



ators £j{0,t) at 



with the eigenvalues aj{t) 



£j{0,t)\aj) = aj{t)\aj). Upon propagat- 



ing through the medium, each pulse experiences a non- 
linear cross-phase modulation. The expectation values 
for the fields are then obtained as 



{£i{z,t)) = Q!l(T)exp 

{£2{z,t)) = a2(T)exp 



1 



27r|a2(r)p 

L6q 
27r|ai(T)p 

LSq 



(14a) 
(14b) 



where ^1,2(2^) — vi^d ± A)LSqz/ {2Tr). These equations 
are similar to those obtained for single-mode jSJ] and 
multimode copropagating fields [l^ . They indicate that 
when the cross-phase modulation is large, upon propa- 
gating through the medium, the phases 



27rsin[61i,2(z)] 



|a2,i(T) 
LSq 



and amplitudes 



ai,2('r) exp <^ -47rsin2[6'i^2(2)/2] 



I«2,i(t)|- 
LSq 



of the quantum fields exhibit periodic collapses and re- 
vivals as 0i^2iz) change from to 27r. In particular, 
when the phase-shift is maximal, 61^2 = 7''/2, the am- 
plitude of the corresponding field is reduced by a fac- 
tor of ri^2 = exp[— 27r|a2,ip/(i<5q)]. On the other hand, 
the maximal dephasing of the multimode coherent field, 

2 = exp[— 47r|a2,ip/(ii5q)], is attained for ^1,2(2) — 
{2n + l)7r (n = 0, 1, 2, . . .), where the phase shift is zero. 

We have thus seen that the behavior of weak quantum 
fields is remarkably different from that of classical fields. 



as in the quantum regime the nonlinear phase shift is 
bounded between zL2TT\a2,i\'^ /{LSq). Only in the limit 
of weak cross-phase modulation 61 2 ^ 1, the quantum 
Eqs. 1141) reproduce the classical result 

(£1,2(2, i)) = ai^2{r)exp[i'q{Ad ± A)|a2,i(T)pz], 

whereby the cross-phase shift grows linearly with the 
propagation distance and can attain large values when 
the field amplitudes are sufRciently high. 

Let us now consider the input state \ipin) = II2), 
consisting of two single photon wavepackets |lj) = 
Y.q^jO'f |0) (j = 1,2). The Fourier amplitudes nor- 
malized as J2q I^JP — 1, define the spatial envelopes 
fj{z) of the two pulses that initially (at t = 0) are lo- 
calized around z — 0, 



l4(z,0)|l,)=^C,V 



qz 



In free space, £j{z, t) — £j(0, r) with t — t — z/ and we 
have (0| £j [z,t)\ Ij) = fj {z — ct) . The state of the system 
at any time can be represented as 



m)) 



Y.^ii{t)\i\)\ii), 

q-q' 



(15) 



e^e2 



from where it is apparent that S^f!^ (0) 

Since for the photon-number states the expectation 
values of the field operators vanish, all the information 
about the state of the system is contained in the intensi- 
ties of the corresponding fields 

{i,iz,t)) = (^i„|£](z,i)4 (z,t) IVin), (16) 

and their "two-photon wavefunction" 0, 



■^,,{z,t;z',t') = {0\£j{z',t%{z,t)\tPin)- 



(17) 



The physical meaning of ^'i-,- is a two-photon detection 
amplitude, through which one can express the second- 
order correlation function G-^^ = 0- The knowl- 
edge of the two-photon wavefunction allows one to calcu- 
late the amplitudes ^'2 of state vector (|15|l via the two 



dimensional Fourier transform of at t = t': 



^ W = l2 // dzdz'^,,(z,z',t)e 



(18) 



We first calculate the expectation values of the inten- 
sities {Ij{z, t)) by substituting the operator solution l|13|) 
into ((TB|l . 



(/,(z,<)) = |/,(-cr)|2 = \f,{zclv, - ct)[^ 



(19) 



where t ~ t—z/vg for Q < z < L. This equation indicates 
that upon entering the medium, as the group velocities 
of the pulses are slowed down to Vg <ti c, their spatial 
envelopes are compressed by a factor of c/vg @. Outside 
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the medium, at z > L and accordingly r = t — L/vg ~ 

{z — L)/c, we have {Ij{z, t)) = \ fj(z + L{c/vg — 1) — ci)p, 
which shows that the propagation velocity and the pulse 
envelopes are restored to their free-space values. 



Consider next the two photon wavefunction 5'^. After 
the interaction, at z, z' > L, we have the general expres- 
sion 



^,j{z,t;z',t') 



U{-CT)f,{-CT'){l + 



Slo 



(20) 



where, as before, t = t — L/vg — [z — L)/c and similarly for r'. For quantum information applications, it makes 
sense to consider the relatively simple case of small magnetic field, such that A, A' ^ A^, where the driving field 
detuning A^ = — W34 satisfies |Ad| < 5uo/2. We thus have 61,2 ~ 6 — rjAdL'^Sq/ {2Tr). Then the equal-time {t — t') 
two-photon wavefunction reads 



*,,(z,z',t) - f,[z + L{c/vg-l)-ct]fj[z' + Lic/i 
f^[z + L{c/vg-l)-ct] . 
f^[z' + L{c/vg-l)-ct] '''' 



-1)- 



ct] 



(21) 



For large enough spatial separation between the two pho- 
tons, such that \z' ^ z\ > Sq^^ and therefore smc[Sq(z' — 
z)/2] ~ 0, Eq. ^ yields 

*y(z,z',i) ~ f,[z+Lic/vg-l)-ct]fj[z'+L{c/vg-l)-ctl 

which indicates that no nonlinear interaction takes place 
between the photons, which emerge from the medium 
unchanged. This is due to the local character of the in- 
teraction described by the sine function. 

Consider now the opposite limit of \z' — z\ <C Sq^^ and 
therefore smc[Sq{z' — z)/2] ~ 1. Then for two narrow- 
band (Fourier limited) pulses with the duration Tp ^ 
\z' — z\/c, one has fj{z)/fj{z') ~ 1, and Eq. l(?T|l results 
in 

^,j{z,z',t) ~ e'''Mz + L{c/vg~l)~ct] 

xf^[z' + L{c/vg-l)-ct]. 

Thus, after the interaction, a pair of single photons ac- 
quires conditional phase shift 6, which can exceed n when 

rsqiy ^ vg\nd\^ 

To see this more clearly, we use Eq. (|18|) to calculate the 
amplitudes of the state vector |'0(i)): 

C!f {t) = e^'^!f iO)exp{^iq + q')[Lic/vg~l)-ct]}. (22) 

At the exit from the medium, at time t ~ L/vg, the 
second exponent in Eq. 1)22(1 can be neglected for all q, q' 
and the state of the system is given by 



When 9 — TT, this transformation corresponds to the 
truth table of the controlled-phase (cphase) logic gate 
between the two photons representing qubits. To- 
gether with the linear single-photon phase shifts (real- 
izing single-qubit rotations), the cphase gate is said to 
be universal in the sense that it can realize arbitrary uni- 
tary transformation [2lj |. 



mL/vg)) = e'^ l^i, 



(23) 



CONCLUSIONS 



In this paper, we have studied a propagation of weak 
probe field through an optically dense medium of coher- 
ently driven four-level atoms in a tripod configuration. 
We have presented a detailed semiclassical as well as 
quantum analysis of the system. One of the conclusions 
that emerged from this study is that optically dense va- 
pors of tripod atoms can support ultrasensitive magneto- 
optical polarization rotation of the probe field and there- 
fore have significant potential for improving the sensitiv- 
ity of Faraday magnetometers to small magnetic fields. 
Another finding is that this system is capable of realiz- 
ing a novel regime of symmetric, extremely efficient non- 
linear interaction of two multimode single-photon pulses, 
whereby the combined state of the system acquires a large 
conditional phase shift that can easily exceed tt. Thus our 
scheme may pave the way to photon-based quantum in- 
formation applications, such as deterministic all-optical 
quantum computation, dense coding and teleportation 
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